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摘 要 : 针对 阅 值 和 阔 值 函数 的 调节 参数 取 值 问题 ， 提 出 一 种 基于 人 工蜂 群 优化 算法 的 带 参 新 冰 值 函数 的 信号 去 唆 算 
法 。 首 先 ， 验 证 带 参 新 闻 值 函数 的 连续 性 、 高 阶 可 微 性 、 参 数 可 调 性 ; 其 次 ， 根 据 最 小 均 方 误差 (MSE) 策略， 利用 
人 工蜂 群 优化 算法 优化 各 分 解 层 的 阅 值 和 调整 参数 ， 得 到 最 优 去 噪 信号 。 最 后 ， 利 用 信 噪 比 《SNR) 、MSE 指标 验证 
信号 的 去 噪 效果 。 实 验 结果 表明 : 人 工蜂 群 优化 算法 选取 的 阅 值 参数 和 新 小 波 阅 值 况 数 可 以 有 效 地 对 带 骂 信号 去 嗓 。 
关键 词 : 小 波 分 析 ; 人 工蜂 群 优化 算法 ; 调整 参数 ; 阔 值 函数 ; 信号 去 噪 
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Signal denoising method based on bee colony algorithm and new threshold function 


Deng Gaofeng, Ye Jincai, Wang Guofu, Zhang Faquan 
(School of Information & Communication, Guilin University of Electronic Technology, Guilin Guangxi 541004, China) 


Abstract: Aiming at the problems of selecting the adjustment parameters in threshold and threshold function, a signal 
denoising method based on artificial bee optimization algorithm and new threshold function is proposed. Firstly, the theoretical 
analysis of the new threshold function is carried out to verify its continuity, high-order differentiability and parameter 
adjustability. Secondly, according to the minimum mean Square error (MSE) strategy, the artificial bee colony optimization 
algorithm is used to optimize the thresholds and adjustment parameters of each decomposition layer, and then the optimal 
denoised signal is obtained. Finally, the simulation experiment is carried out to verify the denoising effect according to the 
signal-to-noise ratio (SNR) and MSE. The experimental result shows that the threshold parameters selected by artificial bee 
optimization algorithm and the new wavelet threshold function can effectively denoise a noisy signal. 
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为 最 基本 的 阀 值 函数 ， 硬 阔 值 和 软 阔 值 函数 都 有 自己 的 缺陷 。 


了 路 


硬 阔 值 函数 的 断 点 问题 使 其 不 具有 连续 性 ， 而 软 阔 值 函数 


到 


品 


as 三 
百 了 次 冰 


号 处 理 技术 的 关键 ， 是 最 大 程度 上 恢复 原始 ”存在 恒定 偏差 的 问题 。 


Pill 


方法 通常 是 基 


是 信 
言 写 的 过 程 ， 它 是 信号 处 理 中 极其 
下 


在 时 频 信号 分 析 中 的 地 位 随 着 信号 的 多 分 辨 分 析 和 局 部 信号 特 。 果 。 但 是 , 对 小 波 系 数量 化 的 空间 还 有 更 进一步 的 提升 。 
征 分 析 的 重要 性 的 提高 而 不 断 提高 。 在 信号 处 理 中 使 用 小 波 变 。 本 论文 提出 一 种 基于 余弦 函数 的 可 导 阔 值 函数 ， 通 过 增 
换 是 空间 域 和 频 域 部 分 变换 的 过 程 。 从 带 噪 的 信号 中 精准 地 提 ” ”因子 来 提高 六 值 函数 功能 的 灵活 性 。 阔 值 的 选取 是 决定 4 
取 有 效 信号 ， 该 过 程 称 为 信号 的 去 品 。 阔 值 函 数 和 浆 值 的 选择 噪 效 果 的 另 一 大 关键 要 素 。 在 小 波 分 解 系数 中 ， 原 始 信号 对 应 


自身 


要 的 一 部 分 。 传 统 的 去 只 为 了 克服 硬 阔 值 和 软 阔 值 函数 的 缺陷 ，2012 年 ，sigmoid 


傅立叶 变换 。 在 过 去 的 十 年 时 间 里 ， 小 波 变换 ”六 值 函数 被 提出 中。 该 六 值 函数 对 小 波 系数 具有 较 好 的 量化 效 


大 


比 ， 


了 一 


[调整 
波 去 


决定 小 波 阔 值 去 噪 效 果 。 的 小 波 系 数 远 远大 于 噪声 系数 中 。 因 此 ， 选 择 合适 的 闵 值 可 将 


闪 值 函数 的 选择 是 决定 小 波 去 噪 效果 的 一 大 关键 要 素 。 小 波 系数 中 的 噪声 系数 有 效 滤 除 。 在 过 去 的 方法 中 ， 闪 值 的 获 
1995 年 ， 美 国学 者 Donoho 提出 的 软 阐 值 和 硬 羡 值 两 种 小 波 立 。 取 主 要 包括 : 基于 邻 域 窗 的 自 适应 阔 值 法 和 经 验 选 取 阔 值 法 ， 


值 去 噪 方法 止 。 这 两 种 常规 阔 值 去 噪 方法 具有 计算 简单 、 有 效 前 者 利用 各 分 解 层 邻 域 小 波 系数 的 相关 性 估计 合理 的 阔 值 欠 


抑 人 


出 噪声 、 保 留 奇异 信息 等 优点 ， 是 一 种 简单 有 效 的 方法 。 作 。 ”后 者 通过 反复 的 实验 ， 最 后 根据 经 验 来 选取 合适 的 闵 值 ， 该 方 
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法 的 工作 量 较 大 。 为 了 克服 人 为 循环 试验 选取 最 优 参 数值 的 缺 
点 ， 并 且 能 够 获得 较 好 的 阔 值 ， 本 文通 过 群 智 能 优化 技术 获取 
合理 的 阐 值 和 调节 参数 。 例 如 粒子 群 优化 算法 、 蚁 群 算法 、 人 
工 鱼 群 算法 等 。 结 合 人 工蜂 群 算法 具有 良好 的 收敛 性 外 、 适 | 
于 多 变量 优化 问题 @、 较 好 的 全 局 搜索 能 力 C8,， 本 论文 采用 人 
工蜂 群 算法 优化 闵 值 和 调节 参数 ， 并 将 该 算法 与 带 参数 闵 值 函 
数 相 结合 ， 


实现 对 带 噪 信号 进行 有 效 的 滤波 。 
1 “小波 立 值 函 数 设计 


小 波 闵 值 包括 硬 闹 值 和 软 闵 值 两 种 ， 这 两 种 方法 的 去 噪 原 
里 都 是 以 特定 闵 值 为 分 界 点 。 不 同 点 是 ， 在 不 同 尺 度 空间 ， 硬 


闵 值 函数 将 绝对 值 小 于 该 闵 值 的 分 解 系数 置 零 ， 同 时 保留 绝对 
值 大 于 该 闷 值 的 分 解 系数 中。 硬 阔 值 函 数 为 
dj (djx|>4 

pa (2-1) 
0 al<4 


而 软 阔 值 函 数 将 绝对 值 小 于 该 冰 值 的 分 解 系数 置 零 ， 大 于 
病 值 的 分 解 系数 通过 相应 算法 滤 除 携带 在 该 系数 的 噪声 数值 。 
软 阔 值 函数 为 


sgn(d jx Wa, 一 4) (dj > 
djx = (2-2) 
0 ldjs| < 
为 了 克服 硬 阔 值 和 软 阔 值 功能 的 缺陷 ， 本 文 提 出 一 种 基于 
余弦 函数 的 新 小 波 阔 值 函数 为 


， (4+1) 
二 一 djr) ld; *A| ， ldjx>=A41 
2 { a Xr+expla/10*(ax|-A) Pal 
0 ， | < 14 
(2-3) 


其 中 dj 为 小 波 系数 ， 儿 为 阐 值 ，Q 为 调整 因子 。 通 


站 整 Q 的 取 值 , 使 该 阅 值 函数 更 好 地 适应 不 同 特性 含 噪 信号 
的 去 噪 。 
很 明显 ， 新 闻 值 函数 由 一 个 指数 函数 构成 。 根 据 基本 函数 
的 特性 可 知 , 新闻 值 函数 同样 具有 连续 性 和 高 阶 可 微 性 。 因此， 
新 闻 值 函数 可 以 克服 不 连续 性 导致 的 伪 吉 布 斯 效应 ， 也 可 以 有 
效 地 抑制 振荡 现象 ， 使 重 构 的 去 噪 信号 更 加 光滑 。 为 了 直接 观 
测 新 闻 值 函数 特性 和 调节 因子 Q 的 作用 ， 简 单 地 取 
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rp 


10，4 = 4。 对 比 不 同 阔 值 函数 图 形 ， 和 


不 同 C 值 的 新 阀 值 函数 图 形 。 如 图 1 所 示 。 图 1 中 ，H-T 表示 
硬 阔 值 函数 ，S-T 表示 软 闵 值 ，N-T 表示 新 闻 值 函数 ，2、1 


20 分 别 表示 C 的 取 值 。 

从 图 1 可 以 清晰 地 看 出 , 当 QQ 越 来 越 小 时 ,新闻 函数 具有 
软 阔 值 函 数值 特征 ， 随 着 C 越 来 越 大 时 ， 新 闻 值 函数 趋向 于 
硬 阔 值 函数 ， 可 知 该 新 阅 值 函数 是 通过 参数 Q 自由 调节 的 软 - 
硬 阔 值 函数 。 
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1 比较 不 同 的 阐 值 函数 和 不 同 参 数 的 新 闵 值 函数 
2 ”人 工蜂 群 算法 (ABC) 


ABC 算法 是 一 种 启发 式 优化 算法 ， 基 于 和 锡 代 的 思想 ， 其 灵 
感 来 自 峰 群 聚 砚 食 行为 00。 标 准 人 造 蜜 蜂 分 为 三 种 , 即 雇佣 蜂 、 
观察 蜂 、 侦 查 蜂 。 雇 佣 蜂 通过 搜索 所 有 的 食物 源 ， 并 对 含有 丰 
富 蜜 量 的 区 域 进 行 搜索 ， 根 据 蜜蜂 之 间 的 行为 交流 ， 食 物 源 被 
观察 蜂 和 侦查 蜂 开 发 并 选 出 最 优 蜜源 。ABC 算法 主要 为 以 下 四 
步 : 初始 化 蜂 群 、 雇 佣 蜂 阶段 、 观 察 蜂 阶段 、 侦 查 蜂 阶段 00。 
2.1 初始 化 蜂 群 

利用 随机 函数 生成 SN 个 nn 维度 的 食物 源 ( 雇 
物 源 产 生 公 式 如 下 : 


Xnin; + rand (0,1)* (x i )) (3-1) 
其 中 i=1,2,3,…， 


蜂 ) ， 食 


ij = 


SN 代表 食物 源 的 个 数 。 


j==1, 2,3,…,n 表示 每 个 食物 源 的 维度 。 Xin ;和 


Xjax,j 分 别 表示 了 维度 空间 的 最 小 值 和 最 大 值 。 


2.2 雇佣 蜂 阶段 


雇佣 蜂 在 各 蜜源 邻 域 进行 搜索 ， 并 且 算 出 每 个 食物 源 的 适 
应 值 与 观察 蜂 共享 。 搜 索 公式 如 下 : 
Vij = Ni, + —)) (3-2) 


其 中 ，9; ;为 [一 1,1] 范 围 的 一 个 随机 数 。r € {1,2…, SN }， 


表示 在 J 维度 中 i 的 邻 域 产生 新 的 食物 源 ， 为 保持 新 食物 源 的 


有 效 性 ，r 去， 通过 贪 禁 选择 , 在 x 与 V 之 间 选 择 好 的 蜜源 ， 


蜜源 的 适应 值 计算 式 如 下 : 


1 
——, F(i)>0 

fitness(i)=41+ F(i) © (3-3) 
1+|F(), F(i)<0 
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其 中 五 (i 为 解 的 


标 函 数值 。 


2.3 观察 蜂 阶段 
观察 蜂 根据 与 密 量 丰富 度 相 关 的 概率 评估 值 选取 最 优 食物 
源 。 这 种 概率 评估 公式 如 下 : 


》 fi(x,) 


其 中 ，fit(x, ) 为 第 i 个 食物 源 的 适应 值 。 蜜 量 丰 富 程 度 与 它 被 


选取 的 概率 成 正比 。 通 过 
源 。 
2.4 侦查 蜂 阶段 


用 轮 盘 选择 方案 选择 概率 估 值 大 的 蜜 


若 在 预定 周期 内 ， 食 物 源 X; ; 没有 得 到 改善 ， 则 放弃 该 食 


物 源 。 侦 查 蜂 变 为 雇佣 蜂 ， 并 且 利 用 公 
域 随 机 产生 一 个 新 的 食物 源 。 


3 ”去 噪 算法 设计 


3.1 目标 函数 
评估 去 噪 后 信号 


\ 式 (3-2) 在 当前 最 大 区 


质量 指标 主要 是 重 构 信 号 与 原始 信号 之 间 


信 噪 比 CSNR ,pw ) 和 均 方 误差 CMSE ) 02， 公 式 如 下 : 
1 这 a 
mse= 志 > (Xx[n]- Xn]] (4-1) 
N 7 一 1 
YY 2 
> XI 
SNR output > 10logio N 和 入 (4-2) 
> x 四 -xD 
1 一 1 
其 中 ， 慎 [四 ] 为 原始 信号 ，X [4] 为 去 品 后 的 信号 。 将 重 构 信 
号 与 原始 信号 之 间 的 MSE 公式 (4-1) 作为 公式 (3-3) 的 目 


标 函 数 ， 构 造 出 本 论文 算法 的 适应 度 函 数 。 
3.2 ”去 噪 参 数 设置 和 算法 流程 

去 噪 算法 分 三 步 

a) 采 用 固定 小 波 变换 〈SWT) 技术 03 对 带 品 信号 进行 小 波 
分 解 ， 得 到 特征 系数 和 细节 系数 ; 
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:ABC 闵 值 量化 re 一 一 闵 值 和 调节 参数 的 优化 


输出 SNR 和 最 小 MSE 


4 
输出 最 优 阔 值 和 调 | 


节 参 数 


输出 去 噪 信号 
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2 去 噪 算法 流程 
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4 ”仿真 实验 


4.1 去 噪 参 数 设 置 

小 波 参数 设置 : 通过 多 次 比较 不 同 滤波 器 和 不 同 分 解 层 数 
的 实验 效果 ， 本 实验 设置 的 分 解 层 为 4， 选 用 Haar 滤波 器 。 

ABC 参数 设置 ， 确定 解 的 维度 等 于 分 解 层 数 与 调节 参数 
个 数 的 和 ， 即 为 D=5， 和 迭代 次 数 为 50， 各 分 解 层 的 阔 值 区 间 为 
[0, 100]， 根 据 图 1 设置 调节 参数 的 区 间 为 [0, 50]， 食 物 源 的 个 
数 SN 为 10， 最 优 解 未 更 新 的 次 数 “limit"* 等 于 0.5*SN*D， 即 为 
25。 
4.2 实验 结果 分 析 

为 验证 提出 算法 的 有 效 性 和 普遍 性 ， 本 文选 
Blocks 基准 信号 作为 测试 信号 ， 分 别 给 这 两 种 基准 信号 添加 不 
同 强 度 的 高 斯 白 噪声 ,产生 SNR 为 10dB、15dB 的 带 噪 信号。 
将 MSE、SNR 作为 去 噪 效果 指标 ,分 别 比 较 硬 (Hard)、 软 (Soft)、 
Sigmoid 阐 值 函数 与 新 闪 值 函数 “New) 的 去 噪 效 果 。 去 噪 数 
据 如 表 1 所 示 。 图 3 为 ABC 算法 收敛 曲线 图 。 图 4 为 SNR 分 
别 为 104B、15dB 时 对 应 的 Bumps 信号 去 噪 图 。 图 5 为 SNR 
分 别 为 104B、15dB 时 对 应 的 Blocks 信号 去 噪 图 。 

表 1 数据 显示 ，Sigmiod 和 新 闻 值 函数 的 去 噪 指 标 均 优 于 
软 、 硬 阔 值 函数 ， 但 是 ， 新 闻 值 函数 可 得 到 更 大 的 信 噪 比 和 更 
小 的 均 方 误差 ， 因 此 ， 对 于 本 论文 的 信号 去 噪 效果 而 言 ， 新 交 
值 函数 的 更 好 。 图 3 的 纵 坐 标 为 MSE， 横 坐标 为 迭代 次 数 。 从 
四 个 子 图 可 以 明显 地 看 出 ， 随 着 迭代 次 数 增加 ， 算 法 不 断 优化 


] Bumps 和 


b) 将 各 分 解 层 阔 值 和 调节 参数 作为 解 的 坐标 分 量 ， 使 用 
ABC 算法 优化 闵 值 和 参数 ， 输 出 最 优 闵 值 和 调节 参数 ; 

利用 各 分 解 层 的 最 优 阔 值 和 调节 参数 进行 信号 去 噪 ， 输 
出 去 噪 后 的 信和 号。 

去 噪 算 法 流程 如 图 2 所 示 。 


冰 值 和 调节 参数 , 从 而 得 到 更 小 的 MSE, 特别 在 算法 的 中 前 期 ， 
收敛 速率 更 快 。 图 4 和 图 5 清晰 地 显示 了 不 同 阔 值 函数 去 噪 后 
的 信号 ， 以 图 4 为 例 ， 软 阀 值 去 噪 后 的 信号 肉眼 效果 最 差 。 硬 
闵 值 函数 去 噪 后 的 信号 依然 存在 大 量 的 噪声 ， 相 比 于 Sigmoid 
浮 数 而 言 ， 新 闻 值 函数 的 整体 去 噪 效果 更 加 的 优秀 。 


本 
可 
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表 1 基于 人 工蜂 群 算法 的 不 同 阔 值 函数 的 去 噪 指标 
有 本 最 小 均 方 误 差 CMSE) 信 噪 比 SNR(dB) 
带 噪 信号 信 噪 比 (d4B) ”原始 信号 
Soft Hard Sigmoid New Soft Hard Sigmoid New 
Bumps 0.2744 0.2413 0.1573 0.1306 16.7401 17.3299 19.1906 20.0016 
10 
Blocks 0.6248 0.5320 0.4305 0.3388 17.5862 18.2167 19.1364 20.1766 
Bumps 0.0894 0.0856 0.0565 0.0516 20.7679 21.7990 23.4495 24.0065 
15 
Blocks 0.3382 0.2499 0.1164 0.1070 20.1838 21.5371 24.4987 25.1829 
1 0.061 一 
09 0.06 1] 
08 0.059 | 
0.058 
D7 有 
0.6 i | 
亚 江 0.056 
拔 05 
时 汉 0.055 
m04 并 0.054 
a 0.053 LL 
0.052 
0.1 0.051 
0 5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 
迭代 次 数 送 代 次 数 
(a)10dB 的 带 噪 Bumps 信号 (b)15dB 的 带 噪 Bumps 信号 
0.9 1.4 
08 ‘2 
则 
S S08 
涉 06 型 
相 相 0.6 
只 os ， 
~ 0.2 
Sa 
3 5 10 15 20 25 30 35 40 45 50 % 5 10 15 20 25 30 35 40 45 50 


选 代 次 数 


(c)10dB 的 带 噪 Blocks 信和 号 


选 代 次 数 


(d)15dB 的 带 品 Blocks 信号 


图 3 人 工蜂 群 算法 优化 收敛 曲线 
原始 信号 带 噪 信号 
- 由 人 八 _ 人员 Eh/ /Na Nome 
~ “本 阅 值 和 和 “ 软 六 值 。 和 人 
o 丰 几 放 只 对 | 从 作 八 人 册 
人 Se 癌 人 ™ Ry 同仁 
对 A 人 A 人 
(a) 
原始 信号 带 噪 信号 
所 10 N 10 
. : 由 八 \ 本 \ 用 . EN/ NN 
200 400 600 800 1000 200 400 600 800 1000 
硬 阔 值 软 闭 值 
jo 10 到 二 10 二 
:人 人 人 贡 对 i 从 作 八仙 
200 400 600 800 1000 200 400 600 800 1000 
Sigmoid 阅 值 , New 浆 值 
a 10 | Nn SS sn 10 人 二 1 
人 六 人员 EA 
9 I pe 400 ee i 1000 I 200 400 、 600 - 而 1000 
(b) 
图 4 Bumps 信号 (a)SNR=10dB (b)SNR=15dB 
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原始 信号 
= 
i Pl 
200 400 600 800 1000 
便 阅 值 
15 > T 
还 10[ | 
| i | ] 
200 400 600 800 
Sigmoid 阔 值 
15 F > T | 
于 位 [ ee = 
到 il 一 N71 1 | | 
200 400 600 800 1000 
原始 信号 
亚 1 r ml 
Te | EE i 
200 400 600 800 1000 
硬 闵 值 
过 1[ | 
EC 一 
200 400 600 800 1000 
Sigmoid 阔 值 
15 T > T 
tf -» = 
EN | . 
200 400 600 800 1000 
图 
5 结束语 
本 论文 提出 了 基于 ABC 算法 的 带 参 新 闻 值 函数 信号 去 噪 
算法 。 通 过 理论 证 明和 仿真 实验 分 析 ， 可 得 出 如 下 结论 : 
a) 本 文 提出 一 种 新 的 带 参 闷 值 函数 用 于 带 品 信号 的 小 波 阐 
值 去 品 。 该 闵 值 函数 在 全 域 具 有 连续 性 、 高 阶 可 微 性 、 参 数 可 
调 性 。 通过 数据 分 析 , 可 知 该 闪 值 函数 可 以 有 效 量化 小 波 系数 ， 


从 而 达到 较 好 的 去 噪 效 果 。 


b) 本 文 将 小 波 阔 值 函数 与 ABC 算法 结合 ， 提 出 一 种 新 颖 


的 去 噪 方法 。 以 阔 值 和 调节 参数 作为 ABC 算法 中 解 的 坐标 ， 


利 ) 


地 得 到 最 优 阔 值 和 调节 参数 。 这 种 方法 上 


j 最 小 化 原始 信号 与 重 构 信号 


Im 


之 间 的 MSE， 可 以 快速 、 有 效 
有 自 适 应 特性 ， 可 以 


克服 传统 阔 值 去 噪 方法 的 不 灵活 缺陷 。 
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